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Abstract 

In this paper we use optimal control to design minimum-time adiabatic-like paths 
for the expansion of a quantum piston. Under realistic experimental constraints, we 
calculate the minimum expansion time and compare it with that obtained from a 
state of the art inverse engineering method. We use this result to rederive the known 
upper bound for the cooling rate of a refrigerator, which provides a quantitative 
description for the unattainability of absolute zero, the third law of thermodynamics. 
We finally point out the relation of the present work to the fast adiabatic-like 
expansion of an accordion optical lattice, a system which can be used to magnify 
the initial quantum state (quantum microscope). 
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1 Introduction 



Quantum science promises computational power unattainable by any classi- 
cal computer [16], as well as unprecedented precision measurement of various 
physical phenomena [23] . At the heart of these important applications lies the 
problem to accurately control and manipulate the states of quantum systems 
For many cases of interest, this control is achieved using adiabatic pro- 
cesses, where the system parameters are changed slowly from the initial to the 
desired final value [12]. If the change is slow enough, the system follows the 
instantaneous eigenvalues and eigenstates of the time-dependent Hamiltonian 
and obtains the desired final state with a good approximation. The inherent 
drawback of adiabatic processes is that they require long times which may 
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render them impractical, since most of the systems are not isolated but are 
exposed to undesirable interactions with the surrounding environment (deco- 
herence) that lead to dissipation [I]. 

Several methods have been proposed to speed up adiabatic quantum dynam- 
ics. Their common characteristic is that they prepare a similar final state as 
the adiabatic process at a given final time (which in principle can be made 
arbitrarily short), without necessarily following the adiabatic path at each 
moment. The corresponding adiabatic-like trajectories are successfully charac- 
terized as shortcuts to adiabaticity. In the method of counter-diabatic control 
[TO] , the applied electromagnetic field restores the adiabatic dynamics of the 
system by suppressing diabatic effects as they are generated. In the method of 
transitionless quantum driving [3], an appropriate auxiliary time-dependent 
interaction is added such that the augmented system arrives in finite time at 
a similar quantum state with the adiabatic trajectory of the unperturbed sys- 
tem. In a closely related method, the trajectory of the shortcut is picked first 
and then the time-dependent interaction generating the corresponding evolu- 
tion is inversely engineered [5] , using the theory of Lewis- Riesenfeld invariants 
|15j . The above methods have been tested experimentally [2"Uf2] and proven 
quite robust to various types of noise [7]. They also share another interesting 
feature: they do not specify a unique shortcut but rather provide entire fam- 
ilies of them [6]. This ample freedom can be exploited using optimal control 
methods to find the shortcuts which minimize relevant physical criteria, like 
time, under realistic experimental constraints [23J. 

In this paper we use optimal control theory to design the shortest adiabatic- 
like shortcut for a potential well with a moving boundary, in the presence 
of restrictions suggested by the experimental setup. This system provides a 
model for a quantum piston and has some interesting applications. Note that 
the shortcuts to adiabaticity for this quantum piston have been studied in [H] , 
as well as their potential applications to the control of many-body quantum 
dynamics. Here we mostly concentrate on the control aspects of the problem, 
but we also highlight some nice applications. In Section [2] we summarize the 
important points of the above work which are essential for the current anal- 
ysis and also formulate the time-optimal control problem for the expansion 
of the piston. We solve this problem in Section [3] and obtain the minimum 
expansion time as a function of the expansion factor of the piston. In Sec- 
tion H] we compare our optimal results with those obtained using the inverse 
engineering method. Using the fact that the minimum expansion time has a 
logarithmic dependence on the expansion factor for large values of the latter, 
we reobtain the cooling rate of a refrigerator [13] , a result which quantifies the 
unattainability of absolute zero implied by the third law of thermodynamics 
|18j . Finally, we discuss how the control problem studied here is related to the 
fast adiabatic-like expansion of an optical lattice with dynamically variable 
spacing [25] . Section concludes the paper. 
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2 Shortcuts to Adiabaticity for a Quantum Piston 



2. 1 Fast Adiabatic-Like Expansion 



Consider the potential 



V p (x,t) 



0, < x < a(t) 
oo, otherwise 



which describes the infinite square well with a moving boundary at x = a(t). 
The evolution of the wavefunction ip(x, t) of a particle trapped in this potential 
is given by the following Schrodinger equation 



2m dx' 



+ V p (x,t) 



(2) 



where m is the particle mass and h is Planck's constant. The wavefunction 
ip is square-integrable on the interval [o,a], with \ip(x,t)\ 2 dx expressing the 
probability to find the particle between x and x + dx at time t. The above 
physical system can serve as a model for a quantum piston. 

We are interested in the expansion of the piston from a(0) = ao > to 
a(T) = ar > ao, at some final time t = T. If the expansion is slow enough 
(adiabatic) then, according to the adiabatic theorem, the solution of ([2]) is 
given approximately by 

oo 

ip{x,t) c n exp [i(f) n (t)]^ n (x,t), (3) 

n=l 



where 



/ 2 • 


U7TX 


/ — TT sm 


[a(t)\ 


ja(t) 



(4) 



are the instantaneous eigenstates of the right hand side in ([2]) [TJ] , 4> n are the 
so-called adiabatic phases, and c n are constant coefficients determined by the 
initial condition. The slower is the expansion, the better is the approximation 
in ([3]). For fast expansion the adiabatic approximation breaks down. 

A method has been proposed recently to accelerate adiabatic quantum dy- 
namics, according to which an auxiliary potential is designed such that the 
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system acquires a similar quantum state with the reference adiabatic path ([3]) 
in arbitrarily short time T [3]. For the case of an expanding piston, it has been 
shown in [5] that if the following potential is added in (121) 



V a (x,t) = h(t)x 2 , 



k(t) = —m 



a{t) 
a(t) 



(5) 
(6) 



then the exact solution of the resulting Schrodinger equation can be expressed 



as 



ip(x,t) = c n exp 



n=l 



exp 



i / E n {t')dt'/h 



(7) 



.ma(t) 2 
l 2ha{t) X 



where 



E n (t) 



n 2 ir 2 h 2 
2ma 2 (t) 



are the instantaneous energy eigenvalues corresponding to the eigenstates 
[14J. The boundary conditions 



a(0) = a , cl(T) = a^, a(0) = a(T) = 0, 



(9) 



ensure the expansion of the piston from a® to in the interval < t < T, 
and that the exact solution ([7j) is similar to the adiabatic one ([3]) at t = and 
t = T. The additional conditions 

fc(0) = k(T) = (10) 



guarantee that the auxiliary potential is active only within < t < T. 



2.2 Inverse Engineering Approach 



The auxiliary potential can be designed using an elegant inverse engineering 
method [5J. The final time T is fixed and the condition (TTOT) is translated 
through to 

a(0) = a(T)=0. (11) 
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A polynomial ansatz satisfying (Q and (TTTj) has been found in [9], a(r)/ao = 
1 + (7 — 1)t 3 (6t 2 — 15r + 10), where r = t/T and 7 = a T /a is the expansion 
factor. Then, the stiffness k(t) of the auxiliary potential can be determined 
from (JO]) 

m 60( 7 -l)r(2r 2 -3r + l) 
1 J T 2 l + (7-l)r3(6r 2 -15r + 10)' 1 J 



where again r = t/T. Note that there is no mathematical limitation on the 
size of T, which can be chosen arbitrarily small in theory. 



2.3 Optimal Control Approach 

In practise, there are always experimental constraints, for example 

-k < k(t) < k , (13) 



which restrict T to some finite value. In such cases, finding the shortest 
adiabatic-like path can be expressed as an optimal control problem. If we 

set 

a(t) frna(t) , , kit) , , 

xi = -^, *2 = J-r-U, u(t) = -f 1 , 14 

a V ^0 a (^J 



and normalize time according to t new = t o \d/T , where T = ym/ko, we obtain 
the following system, equivalent to equation (jSJ) 

ii = x 2 , (15) 

X 2 = — MXi. (16) 

The minimum time adiabatic shortcut, under the constraint f fl3|) . can be found 
by solving the following optimal control problem: 

Problem 1 Findu(t) with —1 < u(t) < 1 such that starting from (xi(0), £2(0)) 
(1,0), £/ie above system with Xi(t) > reaches the final point (xi(T), x 2 (T)) = 
(7,0), 7 > 1, m minimum time T. 

The boundary conditions on Xi and x 2 correspond to those for a and d from 
^2]), while the additional constraint Xi(t) > corresponds to the natural re- 
quirement a(t) > 0. In the following section we solve the above problem on the 
interval < t < T. In order to much the boundary conditions u(0) = u(T) = 0, 
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corresponding to ( ITOj) . the optimal control may be complemented with instan- 
taneous jumps at the initial and final times which do not affect the cost (time), 
see Fig. [TJ in Section HI This approach is similar to that used in our recent 
work [23J2T], as well as in [T3j. Note that the problem of finding the control 
u(t) with < M m i n < u < M max which drives in minimum time the system 
f fT5|) and ffT6]) (the parametric oscillator) from the point (21,22) to the ellipse 
%2 +utx\ = 2Et, with specified oscillator energy Er for stiffness ut, has been 
studied thoroughly in [I]. In Problem [1] the control is allowed to take negative 
values, while the target is a point and not a curve of the state space. 



3 The Shortest Shortcut to Adiabaticity 

The system described by (fToT) and (|T6j) can be expressed in compact form as 
x = f{x) + ug{x), (17) 



where the vector fields are given by 
/ 



22 









, 9 = 









—X\ 



(1* 



and 2 G V = {(21,22) G R 2 : 21 > 0}, u G U = [—1,1]. Admissible con- 
trols are Lebesgue measurable functions that take values in the control set U. 
Given an admissible control u defined over an interval [0,T], the solution 2 
of the system (ll?p corresponding to the control u is called the corresponding 
trajectory and we call the pair (x,u) a controlled trajectory. 

For a constant Ao and a row vector A G (R 2 )* the control Hamiltonian for 
system (fl7|) is defined as H = H(\ , A, 2, u) = A + A[/(2) + ug(x)]. Pon- 
tryagin's Maximum Principle [T7] provides the following necessary optimality 
conditions: 

Theorem 2 (Maximum principle) Let (x*(t), u*(t)) be a time- optimal con- 
trolled trajectory that transfers the initial condition 2(0) = 20 of system (77] ) 
into the terminal state x(T) = xt- Then it is a necessary condition for opti- 
mality that there exists a constant A < and nonzero, absolutely continuous 
row vector function \(t) such that: 

(1) A satisfies the adjoint equation A = —dH/dx. 

(2) For <t <T the functionu h-> H(\o, X(t), x#(t), u) attains its maximum 
over the control set U at u = u*(t). 
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(3) H(X ,X(t),x4t),u,(t)) = 0. 

In the following we use maximum principle to solve Problem [TJ 



Definition 3 We denote the vector fields corresponding to the constant bang 
controls u = — 1 and u = 1 by X = f — g and Y = f + g, respectively, and call 
the corresponding trajectories X- and Y -trajectories. A concatenation of anX- 
trajectory followed by a Y -trajectory is denoted by XY while the concatenation 
in the inverse order is denoted by YX . 

Theorem 4 (Optimal solution) The optimal trajectory for Problem^ has 
the one- switching form XY . The optimal control is 

-1, < t < T x 

u(t) = { X , (19) 

1, T x <t<T x + T Y 



where 



(20) 

TV^i^T). , 21 ) 

The minimum expansion time is 

T = T X +T Y (22) 




Proof. We show first that the optimal control is bang-bang, i.e., alternates 
between the boundary values u = ±1 of the control set. For the system ffT5'|) . 
f lTB]) we have 

H(X , X,x,u) = A + Ais 2 - X2X1U, (23) 



and thus 



Ai=wA 2 , (24) 
A 2 = -Ai. (25) 

Observe that if is a linear function of the bounded control variable u. The 
coefficient at u in H is — A 2 Xi and, since X\ > 0, its sign is determined by $ = 
— A 2 , the so-called switching function. According to the maximum principle, 
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point 2 above, the optimal control is given by u = sign$, if $ 7^ 0. The 
maximum principle provides a priori no information about the control at times 
t when the switching function $ vanishes. Now observe that whenever $(t) = 
— \2(t) = at some time t, then $(£) = —\2{t) = Ai(£) 7^ since the maximum 
principle requires that it is always A = (Ai,A2) 7^ 0. Hence, when = 
it is also 7^ and there is a switch between the control boundary values 
at time t. The optimal trajectory consists of a concatenation of A- and Y- 
trajectories. 

We now move to narrow the candidate sequences for optimality. We show 
first that the concatenation XYX cannot be part of the optimal trajectory. 
Without loss of generality assume that the XY switch takes place at t = 0, 
so A2(0) = 0, while it is also Xi(0) > 0. Along the F-trajectory it is u — 1, so 
from (pM|) and (1231) we find A 2 (t) = — Ai(0)sint. The subsequent YX switch 
should take place at t = 7r, where the switching function $ = — A2 becomes 
zero again. But observe that for u — 1 the state equations (IT5|) and (IToT) 
correspond to a rotation around the origin with period 2tt, so at t = 7r the 
state a; has been rotated by half circle and consequently xi(ir) < 0, which 
is forbidden. We next show similarly that the sequence YXY cannot also be 
part of the optimal trajectory. We assume that the YX switch takes place at 
t — 0, thus A2(0) = 0. Along the A-trajectory it is u — —1, so from (|24p and 
(1251) we obtain \ 2 (t) = — Ai(0)sinht. Observe that $ = -A 2 7^ for t > 0, 
and the subsequent XY switch is not allowed. 

We conclude that the only candidates for optimality left are the XY and YX 
trajectories. Furthermore, it is not hard to see that only the former corresponds 
to expansion (final 7 > 1), while the latter corresponds to compression and 
is rejected. The equations of the A-trajectory starting from (1, 0) and of the 
F-trajectory ending at (7, 0) can be found from (fT5j) and (fTBj) for u = ^fl and 
they are 



xj-x 2 2 = l, (26) 
x? + x^ = 7 2 . (27) 

The switching point satisfies both equations and it is (a/7 2 + 1/ V%, Vl 2 ~ V v2), 
thus using (Tl5l) and (|T6l) the times TV and Ty spent on each segment of the 
optimal trajectory can be easily derived as in f )20|) and f l2Tj) . □ 



4 Results and Discussion 



In Fig. [T] we plot the time-optimal control u(t) from (|T9|) (solid line, T = 
3.4295T ), as well as the control k(t)/k obtained from the inverse engineering 
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Fig. 1. Optimal control (solid line, T = 3.4295To) and inverse engineering control 
(dashed line, T = 6.2511To), under the constraint — 1 < u(t) < 1 and for the 
expansion of the piston by a factor of 7 = 10. Observe that both controls are active 
only within < t < T. 

method (|12p and with minimum duration T under the constraint (1131) (dashed 
line, T = 6.2511To), both for the final expansion factor ax/do = 7 = 10. Ob- 
serve that the control derived from the inverse engineering method is actually 
limited by the lower bound u = —1. 

In Fig. El we depict the corresponding trajectories. Note that the system equa- 
tions (115!) and ( 1T61) can be interpreted as describing the one- dimensional New- 
tonian motion of a unit-mass particle, with position coordinate x% and velocity 
Xi. The acceleration (force) acting on the particle is — ux\. As we can observe 
from Fig. El in both trajectories the particle traverses the same distance x\ 
but along the time-optimal trajectory its speed xi = x<i is always higher. 



In Fig. [3] we plot the expansion time T in units of To = ^Jm/ko as a function of 
the expansion factor 7 for both control policies. This graph provides the speed 
limits of the adiabatic-like expansion under condition (TX3]) - Even if the optimal 
bang-bang controls are not experimentally exactly realizable, knowledge of 
the time-optimal solutions is a useful guide for the design of more realistic 
controls. For example, the abrupt changes of the optimal control shown in Fig. 
[TJcan be approximated by ramps of finite duration [22H3]- ^ n general, more 
complicated constraints on the control or even the state can be incorporated in 
the current formalism using a powerful numerical optimization method based 
on pseudospectral approximations [22J. 

The above results have an interesting thermodynamic application. The quan- 
tum piston using noninteracting particles as the working medium and execut- 
ing the Otto cycle provides a model for a refrigerator, similar to that considered 
in [19]. Following the procedure described in this work, the cooling rate R of 
the refrigerator can be calculated, as the temperature of the cold reservoir r c 
approaches absolute zero. This rate is defined as R = Q/T, the ratio of the 
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Fig. 2. Corresponding trajectories for the control inputs of Fig. [TJ 

heat Q extracted from the cold reservoir on each cycle to the duration T of 
the cycle. We find Q first. When the working medium is in contact with the 
hot reservoir of temperature r^, its internal energy is (ideal gas with one 
degree of freedom). The next step in the Otto cycle is the adiabatic expansion 
of the piston from a to a^. At the end of this process the populations of 
the energy levels are preserved ([3]), while the energies are reduced by a factor 
of 7 2 = a^/a 2 , (jHD) thus the internal energy of the working medium becomes 
r^/27 2 . After the expansion the working medium is brought in contact with 
the cold reservoir and its internal energy is raised to r c /2. The heat extracted 
from the cold reservoir is Q = r c /2 — T] l j2 r y 2 . A necessary condition for op- 
eration of the refrigerator is Q > =>- 7 > \Jrh/T c . For r c — > this implies 
also 7 — > 00, and in this limit the duration of the cycle is dominated by the 
duration of the adiabatic expansion. From fTSOj) and (T2~Tl) we find that for the 
fastest adiabatic-like expansion it is T = Tx +Ty — > In 7 for 7 — » 00, thus the 
cooling rate is restricted as R = Q/T < —t c / lnr c for r c — > 0, and the result of 
[T3] is reobtained. As r c — > the cooling rate approaches zero faster, a man- 
ifestation of the unattainability of the absolute zero (third thermodynamics 
law) PS]. 

We finally mention another possible application of the present work. The op- 
timal control problem solved in the previous section is directly related to the 
fast adiabatic-like expansion of an optical lattice with dynamically variable 
spacing [25]. The connection is that the stiffness u 2 of the external harmonic 
potential, necessary to keep the atoms trapped during the expansion, is related 
to the lattice scale parameter A by the relation u 2 = —A/ A, which is similar to 
(jBJ). As a consequence, the present analysis applies also to this context. Note 
that such accordion lattices are useful since the final lattice spacing can be 
made large enough to be resolved experimentally, for example by imaging of 
the atoms at individual sites. The final quantum state is a scaled-up version 
of the initial state due to the adiabatic-like evolution, thus the optical lattice 
acts like a quantum dynamical microscope [S]. 
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Fig. 3. Expansion time T (in units of To = \fmjko) as a function of the expan- 
sion factor 7, for the optimal (solid line) and the inverse engineering (dashed line) 
strategies. 

5 Conclusion 

In this paper we formulated and solved the problem of minimum-time adiabatic- 
like expansion for a quantum piston, in the presence of experimental con- 
straints. As a result, we obtained the speed limit for this fast quantum driv- 
ing, and used it to rederive an interesting result related to the third law of 
thermodynamics. We also highlighted the possible application of the present 
work to the adiabatic-like expansion of an optical lattice. 
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